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Abstract 

We investigate static spherically symmetric perfect fluid models in Newtonian gravity for baro- 
tropic equations of state that are asymptotically polytropic at low and high pressures. This is done 
by casting the equations into a 3-dimensional regular dynamical system with bounded dependent 
variables. The low and high central pressure limits correspond to two 2-dimensional boundary subsets, 
described by homology invariant equations for exact polytropes. Thus the formulation naturally 
places work about polytropes in a more general context. The introduced framework yields a visual 
aid for obtaining qualitative information about the solution space and is also suitable for numerical 
investigations. Last, but not least, it makes a host of mathematical tools from dynamical systems 
theory available. This allows us to prove a number of theorems about the relationship between the 
equation of state and properties concerning total masses and radii. 


PACS numbers: 04.40.-b, 97.10.-q, 02.90.-bp. 


Keywords: static perfect fluids, stellar models, dynamical systems, 
asymptotically polytropic equations of state. 


‘Electronic address: Mark.HeinzleOap.univle.ac.at 
tElectronic address: Claes.Uggla@kau.se 


1 



1 INTRODUCTION 


2 


1 Introduction 


Spherically symmetric static Newtonian perfect fluid models are the starting point for many discus¬ 
sions about stellar structure and evolution. Among the least complicated models are those described 
by barotropic equations of state p = p(p), where p is the pressure and p the mass density. Within this 
class of models, polytropic equations of state, p = {K, n non-negative constants), have been 

studied thoroughly, and discussions about polytropes can be found in many textbooks on astrophysics 
(e.g., Chandrasekhar Q, Kippenhahn and Weigert ||^, Shapiro and Teukolsky |^). The reason for this 
is at least two-fold: Firstly, they are of physical interest, since they often constitute good approxima¬ 
tions for certain regions of various types of stars. Secondly, the polytropic equations of state give rise 
to symmetries, which makes the equations more tractable. In particular, one can introduce so-called 
homology invariants and thereby reduce the hydrodynamical problem to two coupled autonomous first 
order ordinary differential equations (see e.g., Kippenhahn and Weigert Q, and references therein). 

More realistic stellar models must necessarily go beyond the pure polytropes. The next level of complexity 
is to consider more general classes of barotropic equations of state. In astrophysics numerous equations of 
state describing different kinds of matter, e.g., a stellar plasma or a degenerate ideal Fermi gas, have been 
considered and some of the corresponding stellar models have been studied in detail, although mostly 
numerically. It is worth noticing that most of the equations of state that are of astrophysical interest 
behave like polytropes for low and for high pressures, i.e., they belong to the class of asymptotically 
polytropic equations of state. However, a unified analytical treatment of this important class, comparable 
to that of the polytropes, does not exist so far in the literature, and thus there is a gap that needs to be 
closed. 

In this paper we consider static spherically symmetric Newtonian perfect fluid models with barotropic 
equations of state that are asymptotically polytropic for low and high pressures. This is done by in¬ 
troducing new variables so that a regular set of three autonomous ordinary differential equations with 
bounded dependent variables is obtained. Two of the variables are asymptotically homology invariant. 
Moreover, the variables are chosen so that the equation of state is encoded in a single function appearing 
at only one place in the equation system. For polytropes this function reduces to the polytropic index n. 
Thus the formulation not only naturally incorporates our knowledge about polytropes but also allows us 
to deduce features valid for the general class of asymptotically polytropic equations of state. The main 
goals of the paper are to provide: 


• a framework yielding a visual representation of the solution space, which allows us to obtain qual¬ 
itative information about solutions, 

• a framework suitable for numerical investigations, 

• a formulation that makes a host of mathematical tools from dynamical systems theory available, 
i.e., a framework that makes it possible to prove theorems about e.g., the relationship between the 
equation of state and mass- and radius-properties. 

The outline of the paper is the following: in Section ^ we define and discuss the family of asymptotically 
polytropic equations of state. In Section we introduce new bounded variables, and derive the 3- 
dimensional system of autonomous equations. In Section ^ we discuss the state space associated with 
our formulation and do a local dynamical systems analysis that yields asymptotic expressions for e.g., 
the mass and gravitational potential for the various solutions. In addition we give a theorem describing 
the global dynamics of the state space. This theorem (and most of the following theorems) is proved 
later in Appendix ^ to make the paper more accessible to the general reader. In Section we give a 
number of examples that illustrate the usefulness and pedagogical value of our approach. For instance, 
it is shown that the problem corresponding to Chandrasekhar’s equation of state for white dwarfs has a 
particularly simple and elegant description. We also give a variety of mass-radius diagrams showing how 
extremely complicated mass-radius diagrams can become (at least for unstable models). In Section ^ 
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we state a number of theorems. Firstly, we give criteria for the equation of state that ensure that the 
corresponding models have finite radii and masses. Secondly, we state theorems about how mass-radius 
relationships depend on the equation of state. Most of the proofs are given in Appendix In Section ^ 
we conclude with a discussion about possible generalizations. Appendix ^ contains a brief introduction 
to some relevant background material from dynamical systems theory, and most of the proofs of the 
theorems. Finally, in Appendix ^ we discuss some miscellaneous results. 


2 Equations of state 


Consider a perfect fluid with barotropic equation of state p = p{p). Define 

V = [ dp'p-^{p') . (1) 

Assumption 1. We assume that p is defined for p > 0 and positive for p > 0 and also that rj exists when 

p > 0. 


For static perfect fluids the quantity rj is directly related to the Newtonian potential v. Integrating the 
fluid equation dp + pdv = 0 (Euler’s equation) we obtain rj = vs — v(jp), where vs = n|p=o is the so-called 
surface potential. 

From the above assumption it follows that ry is a monotone function of p and ?7|p=o = 0. Consequently, 
both p and p can be viewed as functions of rj and the equation of state is implicitly given by p(ji) and 
p{r]), where p = p{r])dr]. Moreover, monotonicity of p{p) implies monotonicity of p(ry) and conversely. 

Example. The polytropic equations of state are given by p = where p_ is a positive 

constant and the so-called polytropic index n is a non-negative constant. In terms of p this equation of 
state is given implicitly by p = p-ry” , p = p_7y"+^. 


We now introduce the index-function nijf). Let p(?y) be continuous and sufficiently smooth]^ Now define 


n{r]) 


p dp dp 
p dp ^ dp 


( 2 ) 


Assumption 2. We assume that there exists some Oq > 0 such that n{p) — ng = 0{p°‘°) {p —> 0) and 
some oi > 0 such that n{p) — n\ = 0{p~°‘^) {p —> oo)J^ We also assume that n(p) is on (0, oo) and 
bounded on [0,oo]. 


The above assumption implies that we restrict our attention to asymptotically polytropic equations of 
state, i.e., we consider equations of state that for low pressures (i.e., for p ^ 0; ry ^ 0) take the form 
p{jp) = p_ry"o(l -I- 0{p°'°)), p{p) = jA— _|_ 0[p°'°)), and that for high pressures (i.e., for 

p ^ oo; p ^ oo) take the form p(p) = p+ ?y"i (1 -|- 0{p~‘^^)) , p{p) = p+ry"'^+^(l -I- 0{p~°'^)). This 

corresponds to p(p) OC / Oo + '^) ^ p^l / (ni + l) _|_ Q[jp-ai/(rii + l)'^'^ 

(p ^ oo). 

It is possible to soften the above assumptions, but then it is preferable to incorporate such generalizations 
into the formalism itself. This leads to more cumbersome formulations than the one introduced below, 
and, to show the value of our type of approach, we want to keep things relatively simple in this paper. 
Note that linear equations of state p = const p yield “n = oo” because p does not exist. Equations of 
state that are asymptotically linear are therefore examples that need slightly different treatment than 

^For example, is sufficient. We will not discuss to what extent this restriction can be softened, since we will describe 
below how to handle less restrictive situations like phase transitions. 

^In particular ao (and ai) may be less than 1. Hence, the class of asymptotically polytropic equations of state comprises 
the occasionally considered “essentially polytropic” or “quasi-polytropic” equations of state [iKJ. 
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those presently considered. In Appendix ^ we indicate how to treat equations of state with asymptotic 
behavior p (x p'' > 0 ; where in particular v might take the value 1 ). 

Given n{r])^ an associated equation of state can be written as p(ry) = where the func¬ 

tion 0(77) = 77"“ exp(JJ’ ?7“^ (71(77) — no)drj) (c(0) = 1) and p- = constant. It follows that ^ = 
P-c{r])n{ri)ri^^'^^~^. (The above description is adapted to the low pressure limit, however, if one is so 
inclined, one can equally well adapt to the high pressure limit, and introduce an associated constant p+). 
It follows that a given 77 ( 77 ) gives rise to a I-parameter set of equations of state proportional to each 
other, parametrized by p-. 

Note that 77 ( 77 ) = n = const for the polytropes and that 77 ( 77 ) can be interpreted as a local polytropic 
index. Note also that monotonicity of p{p) is equivalent to non-negativity of 77 ( 77 ). Even though it 
is not necessary, we will for simplicity assume a non-negative 77 ( 77 ) throughout the paper, apart from 
Appendix 


3 Dynamical systems formulation 


In Newtonian gravity the spherically symmetric equilibrium of a self-gravitating perfect fluid is governed 

by 


dm 

dr 


Iirr'^p, 


dp 

dr 



( 3 ) 


where m{r) is the mass inside r; the gravitational constant has been set to 1. Equation (^ is to be 
understood in connection with an equation of state p{p). We will focus on asymptotically polytropic 
equations of state. 


Note incidentally that static stellar models are necessarily spherically symmetric, see for an overview. 
Existence and uniqueness of solutions of (|^ has been proved for rather general equations of state com¬ 
prising the equations of state considered in the present context, see e.g., 0 - 

The above equations (|^) are singular when r = 0, which causes some problems, but even worse is that 
p{p) is not at p = 0 for asymptotically polytropic equations of state. In this section, we will introduce 
new bounded dimensionless variables that lead to a completely regular autonomous system of equations. 
We will do this in two steps since the intermediate step is useful in itself. 


To obtain an autonomous system we ’’elevate” r to an dependent variable and introduce defined by 
^ = Inr, as a new independent variable. We also make a variable transformation from (777, p,r) to the 
following three variables. 


47 rr^p 777 

u = - , q = — , uj = r] , 


( 4 ) 


777 


rr] 


where a is some positive constant (the choice depending on the equation of state, see discussion below). 
Note that u = din 777 /d In r and q = —din 77 /d In r. This yields the following system: 


du , , . dq duj 

— = u{3-u-n{u;)q) , —=q{-l + u + q), —=-aujq. 


( 5 ) 


Here n{uj) is the index-function introduced in (||), i.e., 77(0;) = ^ ^ lr;=wi/“ • 

It is worth noting that for an (exact) polytrope p = 77 ^ (ji a constant), the variable q 

equals q = 77 ^ and hence u and q are so-called homology invariants (see, e.g., Kippenhahn and 
Weigert [||, p.200). In this case the equation for u decouples and leaves a 2-dimensional coupled system 
for u and q. 

Remark. If {m{r),p{r), p{r)) is a solution of (||), then so is (ro 777 (r/ro), ?"(/^p( 7 '/ 7 ’o), r(/^p(r/ro)). However, 
the latter is associated with an equation of state that differs by a factor of from that of the first one, 
i.e., it is associated with {rQ^p{r]),rQ'^p{r])). In the {u,q,Lu;^) formulation, this feature corresponds to 
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translatory invariance of solutions, i.e., the second solution is given by (u(^ — ^o)) — Co))<^(C ~ ^o))) 

where tq = e^°, since ^ = Inr. Accordingly the system (^) describes an equivalence class of proportional 
equations of state {p{r]) = p- c{rf) , p = Jq p{r])dr] \ p- > 0}, parametrized by p_. It is even possible 
to incorporate an additional degree of freedom into the dynamical system by choosing the w-variable 
in a more general way, so that the system describes a two-parameter class of equations of state. 

However, to simplify matters, we are going to take the viewpoint that we always consider a given equation 
of state p{vi) (i.e., we fix p_), and thus that an orbit {('w(^), g(^), a;(^))|^ S R} is associated with a single 
perfect fluid solution (m(r),p(r), p(r)) of (^). 

We now proceed by defining new bounded variables. By our assumptions it follows that oj is positive and 
that u, q are either both positive as well or both negative depending on whether m is positive or negative. 
We now consider the positive mass case, while the negative mass case is treated briefly in Appendix 
The new bounded variables U, Q, H (t/, Q,0 > 0) are defined by: 


U = 


1 + u 


Q = 


1 + 9 


n = 


iO 


1+ UJ 


( 6 ) 


Introducing a new independent variable A according to dX/d^ = (l — U) ^(1 — Q) ^ yields the following 
equations: 


C/(l - [/)[(! - Q)(3 - AU) - n{n) Q{1 - U)] (7a) 

Q(1-Q)[(2C/-1)(1-Q) + Q(1-C/)] (7b) 

-aH(l - H)Q(1 - U). (7c) 

Here, n{Ll) is again the index-function, n{0) = n{ri{Ll)) = ^ ^ |,,(n) = [H/(l —H)]^+ We now include 
the boundaries in our analysis so that our state space consists of the unit cube: [0,1]^. 

To be able to discuss this dynamical system (e.g., to do a fixed point analysis) it is required that n{0) 
is C^-differentiable on [0,1]. This can be achieved by choosing the parameter a sufficiently small. By 
assumption^, n(H) = no+0{fl°‘°^°‘) (H ^ 0) and n(f2) = ni-|-0((l —H)“+“) {ft —s- 1). Therefore, making 
the choice a < min(ao, oi), we obtain dn/(iH|n=o,i = 0 (provided that the 0(-)-terms are well-behaved). 
We are thus mimicking the features of an exact polytrope asymptotically as far as possible. 

Let us now consider composite equations of state where p = p{p) is piecewise on (0,oo). In such 
situations p,m,r are continuous, but dp/dp is not and neither is n(H). It follows from the definitions 
that r], q and hence H, Q are continuous. The mass density p is continuous if p = p{p) is on (0,oo), 
otherwise p makes one or several jumps (a situation that corresponds to phase transitions). The variables 
u and U are continuous if p is, but if p jumps from some value pi to p 2 , then u and U make a jump 
according to ui/u 2 = Ui{l — U 2 )/{U 2 {I — Ui)) = pilp 2 , which follows directly from the definitions of 
u and U. Solutions described by piecewise equations of state on (0, oo) can be obtained by viewing 
the above transformation as a map between state spaces associated with smoothly extended equations 
of state. Note that, in general, if one starts out with a regular solution, then this solution has to be 
matched with a non-regular one associated with the other (extended) equation of state. 

In addition to the dynamical system it is of interest to consider the following auxiliary equations: 


dU 

lx 

dQ 

lx 

dLl 

lx 
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dr 

dm 

lx 


m 


2 


m 


r 


{l-U){l-Q)r 


U{1 — Q)m 
1 / U 


Q 


4-k \ l-UJ \l-QJ 


n 


Ij a 


p-\D) 


1 f U 


Q 


4tt \1-U J \1-Q J 


n 


3/a 


p-^n) 


Q 


i-QJ \i-n 


n 


l/a 


(8a) 

(8b) 

(8c) 

(8d) 

(8e) 


4 Dynamical systems analysis 

An important feature of the dynamical system is that all orbits are strictly monotonically decreasing 
in Lt except ifn = 0, l,(5 = 0or[/ = 1. This fact will be of key importance when we prove the following 
theorem in Appendix A: 

Theorem 4.1. All solutions converge to fixed points when X ±oo except when uq (ni) is equal to 5. In 
this latter case solutions also converge to a 1-parameter set of closed orbits Ci (C 2 ) when X oo(— 00 ). 

Thus the fixed points and the periodic orbits Ci, C 2 describe the asymptotic features of the solutions. It 
is therefore natural to begin with a qualitative analysis of the dynamical system by investigating the fixed 
points. The state space, i.e., the unit cube [0,1]^, and the fixed points of the dynamical system © are 
depicted in Figure |I|. Table p] lists all fixed points together with their local properties, i.e., the eigenvalues 
of the Jacobian of m) at the fixed points. 


Fixed point 

u 

Q 

n 

Eigenvalues 

Restrictions 

Li 

1 

0 

Qq 

1,1,0 


L 2 

3 

4 

0 

Lie 

-Sin 

4 ’ 2 ’ ^ 


A 3 

0 

0 

flo 

3 , - 1 , 0 


La 

1 

1 

Qq 

0 , 0,0 


L 5 

Uo 

1 

0 

1 

1 

1 

T 

0 

no = 0 

Lq 

Uo 

1 

1 

0 , -(1 - Uo) , a(l - Uo) 

ni = 0 

Pi 

0 

1 

2 

0 

nn—3 1 a 

2 ’ 2 ’ 2 


P 2 

0 

1 

0 

^0 5 1 5 ^ 


Pn 

no-3 

2 

0 

1 ^5 - no ± iVb^ , -a(3 

no > 3 

^3 

—2) 

H-no 

Pa 

0 

1 

2 

1 

ni—3 1 a 

2 ’ 2 ’ 2 


P 5 

0 

1 

1 

—ni , —1 , a 


P^ 

ni —3 

2 

1 

1 ^5 - ni ± i'JPj , a/3 

ni > 3 

-^6 

2(ni-2) 

l+ni 


Table 1: Local properties of the fixed points. In the table n(0) = no, n(l) = ni, b = —I — 22ni + Tnf, 
~ (» ~ ~ a; b is negative for Ui < (11 + 8-\/2)/7 and 

positive for ni > (11 + 8^/2)/7. 

Tables ^ and || list from which fixed points (or periodic orbits) the orbits originate and in which fixed 
points (or periodic orbits) they end. If the stable subspace of a hyperbolic fixed point P is rn-dimensional, 
this means that a (rn — l)-parameter family of orbits ends in P. Note that L is a one-parameter set 






















4 DYNAMICAL SYSTEMS ANALYSIS 


7 


Ps 



Figure 1: The compact state space, i.e., the unit cube [0,1]^, of the dynamical system For any index- 
function n{ri) with no > 0 and m > 0 there exist at least 4 fixed points {Pi, P 2 , Pa, P 5 ) and 4 “fixed 
lines”, i.e., lines consisting of fixed points {Li, L 2 , L^, La). In addition, for no > 3 (and/or m > 3), we 
have the fixed point P 3 (and/or Pq). The position of P 3 (Pq) depends on no (ni): P 3 and Pq move along 
the indicated curves as n^ increases. In the picture we have chosen no = 4 and ni = 6 . 
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of fixed points, so that, if the dimension of the stable subspaces of the fixed points on L is m, then L 
attracts a m-dimensional family of orbits. Note that we restricted our attention to limit sets for orbits 
in the interior of the state space in Tables and 


Range of ni 

Fixed point Unstable subspace Dimension 


ni < 3 


P4 < — 4 (ni — 2 )ei -I- 62, 62, eg > 3 

Li < ei, 62 > 2 

L2 <--^n{Ltc)ei+ 2 e 2 + a{Lll-Llc)e 3 > 1 


ni = 3 


P4 — 1-dimensional unstable, 1-dimensional center subspace — 

Li < ei, 62 > 2 

L2 <-^ n(Hc)ei-1- 262 -1-0(0^ - Hc)e 3 > 1 


Pq < 61,62,63 > 3 

Li < 61,62 > 2 

L2 <-^n{Llc)ei + 262 + a{Lll-D.c)e 3 > 1 


3 < ni < 5 



ni = 5 


Pe <61,62, 63 > 3 

Li < 61,62 > 2 

L2 <--^n{Ltc)ei + 262 + nc)e 3 > 1 

C2 -1-parameter family of periodic a-limit orbits- 


5 < ni 


Pe < 63 > 1 

Li < 61,62 > 2 

L2 < —n(Hc)ei -|- 262 + fl(^c ~ ^0)63 > 1 


Table 2: Sources for interior solutions. 


Below we discuss asymptotic physical features associated with the various fixed points and Ci, (72. How¬ 
ever, we already now mention that the solutions with a regular center originate from L 2 . Thus, as is 
well known, there exists a 1-parameter set of regular solutions parametrized by Lie (or equivalently by 
the central density pc for a given equation of state). We refer to this set of orbits as the regular set. 

We then observe that all six faces of the cube are invariant subspaces. Equations (|^) can be solved 
explicitly on the four side faces and we get solutions as shown in Figure ||. 

The = 0 plane (and the H = 1 plane, since it is described by the same equations) display more 
complicated features. Note that the equations for U and Q decouple in equation (^ when one has a 
polytropic equation of state. This leads to a 2-dimensional system which may be identified with the 
Ll = 0 subset (or equivalently the fl = 1 subset). Hence we denote the H = 0 plane as the “polytropic 
subset”. The orbit structure for the polytropic subset is given in Figure |[ 

The polytropic subset is of key importance to subsequent discussions (and highly interesting in itself). It 
is clear that there exist bifurcations when hq = 0,3, 5 (and similarly for rii), which is related to the fact 
that some fixed points become non-hyperbolic. It is worth noticing that the case ng = 0 (corresponding 
to an incompressible fluid) is exactly solvable, as is ng = 5, which is particularly useful. Apart from the 
fixed points, we observe a family of closed orbits Ci centered around P3. The orbits for the ng = 5 case 
can be characterized as the sets {!/, Q} satisfying 

$ = (7(4 - 7Q) - 3(1 - 2 Q) - - Q)^/^ = 0 , (9) 

where the parameter C is required to be (7 > —1/4. This can be shown by observing that ^ = 0. The 
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Range of no 


Fixed point 


Stable subspace 


Dimension 


no = 0 


< C/o(3 — 4C/o)ei + 62 , 63 > 


0 < no < 3 


P2 


< ei, 62,63 > 


no = 3 


P 2 


< 61 , 62,63 > 


3 < no < 5 


Pi < -4(no - 2 )ei + 62 , 63 > 

P2 < 61 , 62,63 > 

P3 < 63 > 


no = 5 


Pi 

P2 

P3 

Cl 


< —4(no — 2 )ei + 62 , 63 > 2 

< 61 , 62,63 > 3 

< 63 > 1 

— 1 -parameter set of periodic w-limit orbits 


5 < no 


Pi < -4(no - 2)ei -I- 62 , 63 > 

P2 < 61 , 62,63 > 

P3 < 61 , 62,63 > 


Table 3: Attractors for interior solutions. 


closed orbits are characterized by C ranging between —1/4 < C < 0; the fixed point P 3 is represented 
by C = —1/4 while (7 = 0 characterizes the solution with regular center, i.e., the orbit connecting L 2 
and Pi- The remaining orbits correspond to positive values of C (compare with Figure ^). It is shown 
in Appendix A that Ci {C 2 ) act as tu-limit (a-limit) sets for orbits coming from the interior of the state 
space. 

Center manifold analysis shows that no interior orbit ends at Pi when no = 3 and that P 4 is a source for 
a 2-parameter family of orbits when ni = 3. For details, see Appendix A. 


To understand the solutions’ mass and radius properties within the state space picture, it is useful to 
consider equations ( | 8 ^ and ( 8 b). Note that they do not involve D, and it thus follows that 2D contour 
plots (corresponding to any section LI = const) yield a picture of where solutions gain radius and mass in 
the present setting (see Figure |[). It is easily seen that no mass is acquired near the side faces U = 0 and 
(5 = 1 , which is in accord with the fact that the fixed points Pi and P 2 are attractors for solutions with 
finite mass (see below for details). In contrast, does not vanish near Pi, but only near P 2 (and P 5 , 
if present), the only attractor for solutions with finite radii. Clearly, orbits spiraling into P 3 must acquire 
both infinite mass and radius. Moreover, as P(A) and Q{X) become asymptotically constant in this case, 
we can specify the asymptotic growth of m(r) as a function of r: m(r) ~ i-(’^o- 3 )/(no-i)^ ggg below. 


Now recall that an orbit stands for a solution (m(r),p(r), p{r)) (r > 0), via (^ and that dlnr = 
(1 — P)(l — Q)dX. The Newtonian potential v{r) is related to r]{r) via v{r) — vs = —? 7 (r), where 
Vs = ' 6 |p=o- For models with finite radius (with surface {p = 0}) vs is the surface potential. At the 
surface the interior solution is joined to an exterior vacuum solution (standard junction). Note that the 
exterior solution is not represented in the interior of the cube. For solutions with infinite radii, p > 0 
everywhere and 65 = 0 . 


We can interpret the different solutions of (i in terms of familiar variables by using the relation between 
T] (i.e, D(A)) and the Newtonian potential and the variable transformations (4). We first analyze the 
orbits in the neighborhood of fixed points from which orbits originate (Table 2), and then turn to the 
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t 

n 



Q-y Pi P2 U -y 


Li Li 


t 

n 




Ps Li 



(a) {/ = 0 


(b) 0 = 0 


(c) C/ = 1 


(d) Q = 1 


Figure 2: Orbits on the side faces, (a) is the U = 0 plane, (b) the Q = 0 plane, (c) the U = l plane and 
(d) the Q = 1 plane. For the latter we had to chose some index-function, the other three are independent 
of n. 


solutions in the neighborhood of fixed points where orbits end (Table ^ . 

Li'. The 2-parameter family of orbits that comes from Li corresponds to solutions with the following 
asymptotic expressions for the potential v(r) and the mass m(r): 

v{r) = Wo + ^ C + b’s + 0{Sr^) (10a) 

m(r) = RC 6r {1 + Sr) + 0{Sr^). (10b) 

Here, vq = ~(yr^)^^“ ^md C > 0 are the two arbitrary constants that characterize the various orbits, 

and Sr = > 0, where R = \j-^C p“^(? 7 )|(-„o) • The orbits that come from Li can be viewed 

as originating from the negative mass cube, as discussed in Appendix A brief investigation of the 
system (^^ and its relatively simple structure (Table ^ shows that there exists a 2-parameter family of 
solutions with a negative mass singularity, for which v{r) = —M/r+Vconst. (where M < 0, Vconst = const) 
when r —!■ 0. For r > 0, these solutions pick up positive mass monotonically. Eventually, when they have 
reached zero mass, they leave the negative mass cube and enter the positive mass cube through Li. This 
happens at the radius r = R (compare with (^). 

L2'. Each fixed point on L2 gives rise to exactly one orbit that moves into the interior of the cube, so 
that the line L2 as a whole generates a 1-parameter family of orbits parametrized by Lie- These orbits 
correspond to the regular solutions of (^, i.e., 

v{r) = Vc +Pct'^ + vs + 0{r'^) (11a) 

Mr) = Pcil + ^nc — r'^) + 0{r'^) (11b) 

3 10 Vc 

when r ^ 0. Here, Vc = )^^° serves as the free parameter, pc = p(77)|(_„^) is the central density, 

and ric is given as ric = n(r])\(^-y^). The regular solutions are certainly the most important perfect fluid 
solutions as they serve as stellar models. 

P4: This point is a source for interior solutions as long as ni < 3. The 2-parameter family of orbits that 
originates from P4 corresponds to solutions asymptotically described by 


= ~~ {^ + Tr^ "1) -b VcouBt+vs + OM) , 

when 

ni 

^2,3, 

(12a) 

M 

v{r) = -(1 - Anp+M r log r) -|- Vconst +vs + 0{r log r), 

r 

when 

ni 

= 2, 

(12b) 

, . const „, , 

v[r) = -^- + Vconst + vs+ 0[-), 

rV— iogr 

when 

ni 

= 3, 

(12c) 
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Figure 4: Subfigure (a) is a contour plot of (b) corresponds to The black regions are values 

of the functions close to one (their maximum value) while the white regions correspond to zero (their 
minimum value). 


in the regime r ^ 0, i.e., the solutions possess a positive mass singularity at r = 0. The constants that 
characterize this 2 -parameter family of solutions when ni ^ 2 ,3 are M > 0 and Vconst ; T is given by 
T = — 47 rp_|_M”i“^(ni — 2)“^(ni — 3)“^. The e appearing in the order-term in ( |l^ ) satisfies e > 0 for 
m < min(2 -|- oi, |) (for oi see Section ||), and otherwise 2 — m < e < 0, which means that 0{r^) 
hides additional singular terms in this case. Like the term T , these terms depend solely on the free 
parameter M. 


For ni = 2, on the space < ei, 62 >, the linearization of (^ has a Jordan normal form and the potential 
looks like in expression (12b). For ni = 3 we have a center manifold and the potential’s leading term 


is as given above. In this case the mass behaves like m(r) 


Y-Iogr 


for r ^ 0. This case represents a 


transition between (1^) and the solutions (p^ generated by Pq for m > 3. 


Pq: There exists a family of orbits that originates from Pq when ni > 3, but one special orbit plays the 
main role. Its potential looks like 


v{r) = (1 -f 0 (r“)) -h vs 


(13) 


as r 


0 . Here, and 


For 3 < ni < 5, (ED is the central orbit of 
In the case (II -I- 8-\/2)/7 < ni < 5, 


477/9+ (ni —1)2 (ni—1) 

a 2-parameter family of orbits asymptotically approaching (|_^ 

we observe a superposition of (H) and dying oscillatory modes. For ni > 5 ( [l^ is the only solution 
generated by Pq. For further details we refer to the treatment of the attractive counterpart of Pq, the 
fixed point P 3 . 


We now turn to the various end points of the orbits: 

L 5 : When no = 0 the 1-parameter set of fixed points L 5 = {(t/p, 1,0)} attracts a 2-parameter family of 
orbits. The corresponding solutions can easily be shown to possess finite radii R and masses M {R and 
M are of course different for the different solutions). To first order the potential has the form 


i{r) = —C 5r + VS , 


(14) 


where Sr = The different orbits are characterized by the two parameters Uq and C > 0. At r = R 



















4 DYNAMICAL SYSTEMS ANALYSIS 


13 


the density is non-vanishing, p = p_, so that solutions of this type can be viewed as asymptotically 
incompressible towards low pressures. The radius R is given by i? = the total mass M 

of a solution is given by M = RC, whereby the surface potential is vs = M/R = C. 

Pi'. Recall that no > 3 is necessary for these solutions to exist. The 1-parameter family of orbits that 
converge to Pi corresponds to solutions with infinite radii and possess the following potential, mass, and 
density as r ^ 00 : 


(l-^-^C^r3-"'>) + 0(.) 

(15a) 

r no — 6 

M(1 r^-^°) + 0{-) 

(15b) 

no-3 

p_M"0r-"“(l-hO(r3-’"«)). 

(15c) 


Here, C denotes C = ^ and M > 0. The solutions (15) thus have finite masses M (different 


for different solutions) although they have infinite radii. For the special case of an exact polytrope with 
n{r]) = 5 the associated regular solutions are given by v{r) = —M/^J ^ p-M^ + r^. 


P 2 : When riQ > 0, the fixed point P 2 acts as an attractor for a 2-parameter family of orbits. Since these 
orbits correspond to solutions with finite masses and radii, they play the most prominent role in our 
considerations. 


v(r) =—C Sr — C Sr^ + vs + 0(Sr^^°‘) (16a) 

/ rjno IJ-l \ 

m(r) =RC i 1--f O(-) (16b) 

V no +1 ) 

p(r) = -hO(-). (16c) 

As before, Sr stands for Sr = The constants C > 0, P > 0 are arbitrary and a = min{l,no}. 

Combining the earlier asymptotic expression for p(r]) with the above equations, one can show that the 
radii and masses of the solutions that end at P 2 are given by R? = 4 ^^ ^ = 4 ^^ ^. 

P 3 : For 3 < no < 5, the fixed point P 3 attracts a single orbit that corresponds to a solution with infinite 
mass and radius. For r —> 00 we have 

v(r) = 0(r-“)) (17a) 

m(r) = (1 0(r-“)) , (17b) 

no - 1 

where a = and AT = 4 ^^ For np > 5 the solution (|^) becomes the center of a 2-parameter 

family of orbits that oscillate around ( 0 ) and converge to it when r —> 00 . 

v(r) = (1 -I- Hr~°'^ osc(a 2 logr) -|- 0 {-)) , (18a) 

where osc(x) = (^102 — ^ 201 ) cosx -I- (fcioi -I- k 2 a 2 ) sinx 

m(r) = —^— Kr^^o-^yCo-i) (j^ _|_ bsc(a 2 logr) -|- O(-)) , (18b) 

no - 1 

where osc(a;) = (1 — ai) osc(a;) -I- 02 osc'(a:) 

In (|^) K is again defined as in (|^; furthermore, ai = 2 (n^-i) > «2 = 2 (^- 1 } , and iJ = 16 
The constants fci and ^2 parameterize the various solutions. Note that v(r) and m{r) are only of the 
form (p^, if the condition 4ao > (no — 5) is satisfied (whereby a > Oi); else the oscillations are dying 
out more rapidly than the 0{r~°‘) term in ( 0 - 
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When no = 5 (ni = 5) the periodic orbits Ci (C 2 ) constitute limit sets for interior orbits and it is possible 
to obtain asymptotic expressions. We refrain from giving such expressions, but to obtain them one uses 
the known exact solutions for polytropes with n = 5, whose asymptotic expressions directly yield the 
desired expressions for more general equations of state since dn{Ll)/dPl — 0 when ^ 0. 

Remark. We ignored those fixed points which are not mentioned in the list since we are only interested 
in interior orbits. Only L 4 requires a closer investigation, since the linearization of ( 0 ) at L 4 is zero. 
However, considerations in Appendix A show that no interior orbits converge to L 4 for A —> ± 00 . 

Remark. Finally note that the orbits on the side face U = 0 (see Figure ||) correspond to the vacuum 
solutions v{r) = —^ + C + vs {M > 0), where (7 = 0 corresponds to the orbit P4 ^ Pi, while (7 < 0 
(C > 0) characterizes the orbits on the left (right) hand side. The orbits on the side face Q = 0 correspond 
to solutions with constant potential. The orbits on U = 1 and Q = 1 are unphysical since dlnr = 0 in 
these cases. 


5 Examples 

In this section we give some examples and discuss some features concerning the relationship between 
the equation of state and various models’ radii and masses. We will come back to this issue in the next 
section where we give some general theorems. 


5.1 Mass-radius features 


Recall from the discussion of the fixed point P 2 that it acts as an attractor for a 2-parameter family of 
orbits and thereby generates a 2-parameter family of finite perfect fluid solutions (cf. dH))- An orbit of 
this family can be uniquely characterized by two parameters (7 > 0 and D > 0, where C = limA—.cx 3 

and D — lim;^^oo . The associated physical variables display the asymptotic behavior (H) asr ^ R, 
and for the finite radii R and total masses M the relations 


R^ 


1 C 
47rp_ D 




1 (73 

47rp_ D 


(19) 


hold. We may thus refer to ( |l^ as a universal (M, i?)-relation. 

Consider a 1-parameter family of solutions, described by a parameter s (e.g., the regular solutions with 
s = Pc the central density). We then have (7(s), D{s), and hence R{s) = R{C{s), D{s)); M{s) = 
M{C{s), D{.s)), whereby we obtain a curve in a (M, i?)-diagram for the particular family under consid¬ 
eration. This can be visualized as making a particular 1-dimensional cut through the graphs R{C, D) 
and M{C, D), defined by (H). 

Equation (|l^ can be brought into a more illustrative form by introducing cylindrical coordinates around 
P 2 , i.e., set U = esincj), 1 — Q = ecoscf), Lt = h‘^ (with e small). Then the asymptotic expansion around 
P 2 yields. 


r 


2 


47rp_ 


tan())-h(7(e“), 


—-tani^ (1-|-tan^ ()))-I-(7(e“) 

47rp_ e‘‘ 


( 20 ) 


as A ^ 00 , and, accordingly, 

M = (47rp_)-i/(”«-i) lim [e-i(tan((.)i/("“-i)(cos<^)-i] . (21) 


Equation ( 0 ) and related formulas are useful in e.g., numerical computations: take a certain small fixed 
e and follow an orbit attracted by P 2 to its intersection with the small cylinder. Then the values of 4> 
and h of the point of intersection yield M and R. 
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Note that (^) and (H) automatically yield the standard relations (see, e.g., |^) for the polytropes 
n(rj) = n: Since (Q) decouples for polytropic equations of state, all corresponding regular orbits intersect 
the small cylinder at the same value of </>. Moreover, dzi) can be integrated, and expressed in rj we obtain 
r?(A) = rjcexp ( — Qp(A)(l — Up(A))dA ), where (C/p, Qp) is the solution of (f^ ) and At some 

value A = A the orbit intersects the cylinder, where A is independent of r/c, whence we conclude that 
h oc 77c. Therefore, for the polytropes, 

= const = const r]^~^ M = const (22) 

where the constants only depend on p_ and on n (but not on pc). {M, i?)-diagrams for some polytropes 
are displayed in Figure |^. 




(a) n = 0.25 (b) n = 1 (c) n = 1.5 



(d) n = 3 


(e) n = 4 


Figure 5: The mass-radius diagram for polytropic equations of state with different polytropic indices n. 
We use bounded variables, M = (M/M)(l -|- M/M)~^ and TZ = {R/R){1 + R/R)~^, where M and R 
are typical values. 


5.2 Chandrasekhar’s equation of state 

The equation of state of a completely degenerate, ideal Fermi gas reads (see e.g., 

p{x) = CpX^ , p{x) = Cp[x\/ l + x"^ (2x^/3 - 1) -I- log(a; + \/1 + x"^)) , 


(23) 
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where Cp and Cp are constants. The equation of state is given in implicit form and the parameter x is the 
“dimensionless Fermi momentum”. In the non-relativistic limit, x the pressure p{x) oc x^ , and in 
the ultra-relativistic limit, x 1, p{x) oc x^. Hence the equation of state behaves like a polytrope 
with no = 3/2 (ni = 3) for low (large) densities. 

In terms of rj the equation of state and the polytropic index function are represented by simple 
expressions. For convenience, we use a dimensionless variable t], which is just a rescaled version of our 
original 77 -variable according to 77 /cp, where Cp = {8/3){cp/cp). We obtain 


piv) =p.{v{l + ^)f\ (24) 

and 

^ 

This latter expression makes it evident that we deal with an asymptotically polytropic equation of state 
for which no = 3/2, ni = 3, and that oo = oi = 1. 


Clearly, as is well known, all stellar models with a completely degenerate ideal Fermi gas have finite radii 
and masses since ng = 3/2 (compare with the theorems in the next section). 


The compact state space (U, Q, fl) provides a picture of the solution space of a given equation of state. 
In Figure projections onto the If = 0 and U = 0 planes of the regular “Chandrasekhar solutions” are 
shown. In Figure |6(a) , for small central densities the orbits coincide with the n = 3/2 polytropic orbits. 
With increased central density the orbits bend down until, in the limit of infinite densities, the orbits 
approach the n = 3 polytropic solution. 


P 2 


I/4 L3 Pi P5 




(a) Projection to ft = 0 


(b) Projection to f/ = 0 


Figure 6: Solutions associated with Chandrasekhar’s equation of state projected down to = 0. The 
topmost orbit coincides with the polytropic orbit for n = 3/2, the lowermost orbit coincides with the 
polytropic orbit for n = 3. 

Also the corresponding mass-radius diagram “interpolates” between the polytropic ends. For small 
central densities we observe that M and R behave like for a polytrope with index n = 3/2, while for high 
central densities the mass approaches the Chandrasekhar limit and thereby becomes independent of Pc, 
analogous to a polytrope with index n = 3. For the mass-radius diagram see Figure |^. 


5.3 Composite equations of state 

Since the equation of state only enters the system of equations through the function n(H), one can view 
the solutions as a ’’changing” polytrope. Figure ^ gives an intuitive idea of how solutions behave in the 
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Figure 7: The mass-radius diagram for Chandrasekhar’s equation of state. We use bounded variables, 
Ai = (M/M)(l + M/M)~^ and 7Z = (P./R)(l + R/R)~^, where M and R are typical values. 


compact state space. It is possible to approximate an equation of state by discretizing n(fl) and thus to 
view it as a composite equation of state. 

We will now give some mass-radius diagrams for simple composite equations of state: consider equa¬ 
tions of state that correspond to functions n(fl) described by two constant values, no,ni, and a jump 
from no to ni at a given LI = Clj. In this present case U,Q,d are continuous (compare with the remark 
at the end of Section ||) and the map between the state spaces associated with the polytropes no,ni 
can be viewed as gluing the state spaces together at IIj. Hence a regular solution that corresponds to a 
composite equation of state follows the polytropic flow associated with the index ni until Cl has reached 
Lly At this point the polytropic flow associated with uq takes over. 

Figures || are projections to the = 0 plane. Regular solutions proceed along the spiral (ni-flow) up to 
some point, at which they switch to the np-flow. If this point lies on the curve that ends at Pi, then the 
solution has infinite radius but finite mass (see (H))- 

We observe that when no is chosen between 3 < no < 5 and m > 5, a subset of the following possible 
scenarios occur: 

• The two spirals do not intersect and thus every regular solution has a finite mass and radius. See, 
e.g.. Figure 8 (a). 

• Finitely many solutions with infinite radius - but with finite mass - are embedded into a continuum 
of finite solutions. Compare with, e.g.. Figure 8 (b). 

• A continuum of finite solutions contains one solution with infinite radius and mass. This solution 
represents the accumulation point for an infinite number of solutions with infinite radii and finite 
masses. This happens if the (projection of the) ni-spiral passes through P 3 . See, e.g.. Figure 8 (c). 

• Infinitely many solutions with infinite radii and finite masses are embedded into a continuum of 
finite solutions, but without an accumulation. This is the case if Pq possesses {U, Q) coordinates 
lying on the no-spiral. 

• For a certain special choice of no,ni the last two scenarios can occur simultaneously. 

Using relations describing the polytropic spirals like (|^), it is possible to give a qualitative answer to 
the question which central values yield infinite solutions. 
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(a) no = 4, ni = 7 


(b) no = 4, ni = 6.2 


(c) no = 4, ni = 5.88 


Figure 8: Composite equation of state 71(77) = no for rj < rjj and 71(77) = for rj > rjj. For (a), (b), and 
(c), tiq = 4 while ni varies. The orbits are projected down to fl = 0. A regular solution starts from L 2 
and proceeds to some point (depending on rjc) on the rii-orbit (where fl (ji) has reached flj (77j)). Here 
the orbit is matched with the rio-flow, which is subsequently followed. A regular solution with infinite 
radius and finite mass is constructed by gluing the regular orbit, which originates from L 2 , to the orbit 
that ends at Pi. 



Figure 9: Composite equation of state 71 ( 77 ) = ''^0 for V Si Vi s-nd n{r]) = 711 for ij > rjj. tiq = 6 for (a) and 
(b), Til varies. Orbits projected down to fl = 0. Compare with Figure 0. 
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For composite equations of state with ng > 5 the following is possible: 

• Every regular solution has an infinite radius and mass (Figure 9(a)). 

• There exists one solution that has infinite radius, but finite mass. This solution separates a con¬ 
tinuum of solutions with infinite radii and masses from a continuum of solutions with finite radii 
and masses. See, e.g., Figure 9(b). 

As seen from Figures ^ and |lT|, the mass-radius diagrams can become quite complicated already for 
this simple class of equations of state. However, by juxtaposing Figure and Figure § the qualitative 
features of the mass-radius diagrams can easily be explained. 

Note that it is easy to construct smooth equations of state displaying similar features. One could, e.g., 
study equations of state of the Chandrasekhar type ( p^ with general constants. One can also choose to 
use some suitable regularization at the jump to obtain solutions, with a smooth equation of state, that 
display similar features. 


6 Qualitative results 


We start this section by deriving a number of inequalities that give information about the solution space 
for various equations of state. 

For regular solutions it is clear that p{r) < p{r), where the average density p is defined by p = 3m/(47rr^), 
if dp/dp > 0 (or equivalently n(H) >0). In our approach this inequality translates to the simple statement 
that U < 3/4. This follows from the fact that u = 3p/p. An alternative proof is obtained by first noting 
that dU/dX = —3n(H)Q/64 on the U = 3/4 surface, and that dU/dX < 0 if 17 > 3/4 and n{fl) > 0. We 
then note that the eigenvector for the transversally hyperbolic set L 2 has the components (—3n(Hc)/40, 2) 
in the (C/, Qj-plane. Hence regular orbits start on the inside of the surface or on it if n(H) = 0; in the 
latter case they continue on the surface as long as n(H) = 0. Thus, because of the inequality, it follows 
that they cannot pass into the region where U > 3/4. 

Another inequality follows from that Pc ^ p > p, when dp/dp > 0 (n(H) > 0), and the definition of U. 
It is given by p{Vt)/pc < C//(3(l — U)). Combining this with equation (||) yields 


U > 


Q 


47r \1 — Q J \1 — H 




l/a 


.,-1 


m > — 


Q 


In \1-QJ Vl-H 


n 


3/a 


(26a) 

(26b) 


As another example, which will be used in proofs of theorems below, consider a polytrope, characterized 
by n, and another equation of state, characterized by the index function n(Q), such that n > n{Q) (or 
n < n{n)) for 0 < H < Hmax- Now consider a 1-parameter set of solutions of the polytrope, e.g, the 
regular solution set for which 0 < He < Hmax- Since the equation for H decouples in the polytropic 
case, it follows that the polytropic invariant set is described by a particular surface (17,(5) = 0 in the 
state space U,Q,fl, and that dlf/dX = {dl>/dU){dU/dX) {dl!/dQ){dQ/dX) = 0 on the surface. We 
restrict ourselves to subsets that satisfy dU/dX < 0 and dQ/dX > 0 in the interior state space. The 
regular subsets for polytropes with n < 5 are particularly interesting subsets of this type (this will be 
shown below; an intuitive picture of what types of sets that satisfy this requirement can be obtained by 
considering Figure ||). From the above it follows that {d'^/dU) and {d'^/dQ) have the same (non-zero) 
sign. We can always choose so that the sign of (i94>/i9{7) is positive. 

Let us now consider the described ’’polytropic surface” 'E{U,Q,fl) = 'E{U,Q) = 0 in the state space of 
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Vc 

(a) no = 4, ni = 6.2 



Vc 

(c) no = 4, ni = 5.88 



Vc 

(b) no = 4, ni = 6.2 



(d) no = 4, ni = 5.88 




(e) no = 4, ni = 5.88 (f) no = 4, ni = 5.88 


Figure 10: Composite equation of state 71(77) = no for rj < rij and 77(77) = ni for 77 > Tyj. tiq = 4 and 
77i = 6.2 or 77i « 5.88. The diagrams display the radius R and the total mass M of regular solutions 
as functions of the central value rjc- For tzi = 6.2, R becomes infinite for two values of rjc, while M is 
finite, although the mass can become arbitrarily large. For tzi « 5.88, we observe an infinite number 
of solutions with infinite radii and finite masses. These solutions accumulate at a certain value of r]c, 
where both R and M are infinite. By comparison with Figure 8(b) and 8(c) the qualitative features 
of the R- and M-diagrams can easily be explained. In the diagrams bounded quantities are used: 
fjc = (77c/10)(l + 77c/10)“^, M = log(l + M)(l + log(l + M))“\ and TZ = log(l + R){1 + log(l + R))~^- 
Note that TZ = 1 corresponds to an infinite radius. 
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M 


M 


(a) no = 4, ni = 6.2 


(b) no = 4, ni = 5.88 


Figure 11: Mass-radius diagrams for composite equations of state, i.e., n(ji) = no for rj < rjj and n(ry) = ni 
for r] > rjy For no = 4, ni = 6.2 the diagram can be understood as follows: For small rjc (M, R) is located 
in the right upper corner of the diagram. With increasing a point {k{M), k{R)) « (0.35,1) is reached, 
then there is a jump to a small loop graph in the right upper corner. Having passed through the loop 
there is a jump to a point {k{M), k{R)) « (0.45,1), and subsequently the curve goes down until it ends 
up in a (M, i?)-spiral (not completely represented in the diagram). For ni « 5.88 the (M, i?)-diagram 
becomes more complicated. As is only indicated, in the upper right corner an infinite number of small 
graphs accumulate. In the diagrams bounded quantities are used: A4 — log(l -I- M)(l -|- log(l -I- M))~^ 
and TZ = log(l -I- i?)(l -I- log(l -b R))~^- 


the other equation of state. An evaluation of the derivative of 4* on the surface yields 
\ / n(n) n{n) \d^Jn 



= (|^)c/(l-C/)"Q(n-n(H)), 


where we have made use of a pointwise comparison of the two flows and where the one associated with 
n is zero on the surface. The above quantity is positive (negative if n < n(Ll)) in the interior state space 
and acts as a ” semi-permeable membrane” for solutions with 0 < < flmax- Note that if = Ij this 

surface will provide a restriction in the state space for all solutions. As an example, we now consider 
n = 5 with Hmax = 1, which plays a key role in several of the theorems concerning models with finite 
radii below. 


As mentioned earlier (see the regular subset for n = 5 is characterized by 


3(1-2Q) 
4-7Q 


(28) 


Now define 

<^(U,Q,fl) = (4-7Q)U-3(1-2Q) . (29) 

We now consider the surface 4)(U,Q,fi) = 0 inside the state space, and note that the surface has an 
outward pointing normal. The flow through the surface is given by the derivative of <i> on the surface. 
Equation (^ yields 




surface 


3(1 - 2Q)(1 - Uf Q(5 - n(fi)) . 


(30) 
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Note that Q < 1/2 on the surface, and hence that the sign in the interior state space is determined by 
5 — n(f2), in accordance with When n{Ll) < 5 (> 5) the derivative is positive (negative) and the 
vector field is pointing ’’outwards” (’’inwards”), as can be intuitively deduced by considering Figure ^ 
when one superimposes the regular orbit for u = 5 onto the figures for the other polytropic indices. We 
also note that by continuity the equilibrium point P 3 (Pq) is inside the intersection of the surface with 
the n = 0 (n = 1) plane. Thus for example, it follows from this inequality, together with the eigenvector 
properties associated with L 2 , that the regular surface of n{Ll) < 5 (> 5) is ’’outside” (’’inside”) the 
regular ”n = 5-surface” (once orbits are outside (inside), which they are in the neighborhood of L 2 , they 
are prevented by the inequality to pass through the surface). 

Moreover, since it follows from equation (||) that dU/dX < 0, dQ/dX > 0, if n(fl) < 5 , outside the regular 
n = 5-surface, the above ’’outside” statement is also true for u if 5 > n > n{Il). 

Let us now list a number of theorems, relating equations of state to properties concerning total masses 
and radii of various solutions (primarily the regular ones). For the proofs (which are mostly given in 
Appendix ^ we can make use of methods provided by dynamical systems theory. This is a considerable 
advantage compared to the cumbersomeness of several proofs in the literature; for theorems of a similar 
type as the theorems below see, e.g., §, [Q, [Q, |^, ||ll|. In addition to our unified approach to 
already known statements we will also prove new results. 

Theorem 6.1. (Finiteness of perfect fluid solutions). Consider an asymptotically polytropic equation 
of state p{p) = iyp"o/("o+i) (1 -I- (p 0^ with asymptotic index no < 3. Then any 

corresponding perfect fluid solution, either regular or non-regular, has finite mass and radius. 


Proof. We distinguish two cases, no < 3 and np = 3. For no < 3 the proof is trivial. According to 
Theorem O only fixed points can be attractors on 12 = 0 in this case. Combining this with Table H, it 
follows that there exists a unique attractor for np < 3, namely L 5 for np = 0 and P 2 for 0 < np < 3. 
These fixed points correspond to solutions with finite radii and masses (compare with ( 0 ) and (|l 6 |)). 

In the case np = 3 the argument must necessarily also involve Pi, which is no longer hyperbolic, but 
possesses a I-dimensional center subspace. Therefore it cannot be excluded that Pi is an endpoint of a 


family of orbits that comes from the interior of the state space. However, in Appendix A.2 it is shown 
that this is not the case. Consequently, it follows that all orbits end at P 2 , and hence Theorem | 6 .l| is 
established. ■ 


Theorem 6.2. Consider an asymptotically polytropic equation of state pfq) with index-function nirf) < 5 
V 77 , n{r]) < 5 for arbitrarily small p. Then the regular perfect fluid solutions, and the solutions that possess 
a negative mass singularity, have finite masses and radii. 


Proof. Equation (g^ defines the surface of polytropic regular solution of index 5. It follows from the 
properties of the eigenvectors of the transversally hyperbolic line P 2 and from (|3^) that this surface 
shields the regular solutions, for any n(f2) < 5, from the fixed point P 3 (when it exists). Taken together 
with Theorem 4.1 and the properties of Pi, it follows that P 2 if np ^ 0, and P 5 , if np = 0, are the only 
limit sets when A —*■ 00 for the regular solutions, and thus that they have finite radii and masses. ■ 


Theorem 6.3. For any np > 3 there exist asymptotically polytropic equations of state (having np as 
their asymptotic index), such that some of the regular perfect fluid solutions possess infinite radii. 


Proof. When ni = 5 -I- e, for some small e, the regular subset for ni intersects the surface associated with 
orbits that end at Pi for an equation of state with polytropic index np. Hence one can obtain an example 
displaying the desired feature by constructing a composite equatio n of state. For further discussion, and 
a picture of solutions exhibiting this type of property, see Section ^ and Figure |[ ■ 


Hence Theorem 6.1 is sharp. From Theorem |6.2| we know that an equation of state that yields regular 
solutions with infinite radii has to satisfy n(f2) > 5 in some range for Q. However, this condition is not 
sufficient, since there exists some equations of state that satisfy n{Lt) > 5 in some range for H, for which 
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all regular solutions have finite radii. We refer to the examples given above for composite equations of 
state. 

The dynamical systems formulation is also a powerful tool to investigate the qualitative properties of the 
mass-radius diagram. A spiral structure of the (M, ii)-diagram has been observed in several situations 
and its existence has been proved rigorously in some cases (see, e.g., [^). The following theorem 
formulates a criterion when such a structure occurs and gives information about its concrete form. 

Theorem 6.4. (Spiral structure of the (M, R)-diagram). Consider an asymptotically polytropic equation 
of state with asymptotic indices 0 < no < 3 and ni. Then the mass-radius relation for high central 
pressures possesses a spiral structure, if and only if ni >5 with the spiral given by 


( RiVc) 

V Mip,) 


Ro 

Mo 


+ 


1 1 

— ) e J 7(72 log - ) 6 
VcJ V 



(31) 


where Rq and Mq are constants, B is a non-singular matrix, and b a non-zero vector. The matrix 
J{gS) G SO(2) describes a rotation by an angle ip, and the constants 71 and 72 are given by 


71 = (m - 5) i 72 = i ^6 . 


(32) 


Recall from Table that \/b = 

Proof. See Appendix ■ 

Theorem 6.5. (Polytropic behavior of the {M, R)-diagram). Consider an asymptotically polytropic 
equation of state with asymptotic indices 0 < no, and 0 < ni < 5. To first order, the mass-radius 
relation for high central pressures is approximated by the mass-radius relation for an exact polytrope with 
polytropic index ni, i.e.. 


RiVc) = Rp-.m iVc) C[n{v)] M(ric) = Mp,n^ {rjc) C[n{ri)] 


(33) 


where Rp-.miVc) and Mp-mipc) are the radius and the mass for an exact polytrope with index ni (compare 
with i.e., 

Rp-.miVc) = Cfi(ni) Mp,„fir]c) = cm(?^i) , (34) 

and is a functional acting on the index-function 71 ( 77 ), by 

1 


C[n{r])] = exp I i J 


njv) - no 1 f n{r]) - 

V 2 7 77 

0 1 


dr] 


(35) 


Proof. See Appendix 


7 Concluding remarks 

In this paper we studied static spherically symmetric perfect fluid models in Newtonian gravity with 
asymptotically polytropic barotropic equations of state. We developed a dynamical systems formulation, 
which yielded both a pictorial representation of the solution spaces, and a tool making it possible to 
derive theorems. 

The framework was also used for numerical purposes, and exhibited certain advantages, in particular 
in connection with bifurcation phenomena. Such bifurcations cause numerical problems, which can be 
circumvented by use of approximate solutions near fixed points. 
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In the paper, the approach is to some extent tailored to the family of equations of state that are considered. 
However, the underlying principles have a broader range of applicability - in the present context of 
Newtonian perfect fluids, but also in other areas. 


Let us first comment on the connection between the choice of variables and the asymptotic features of 
the equation of state. The present choice of the variable a; = 77“, in (^, is adapted to the asymptotic 
features, 71(77) = + 0(77““), 71(77) = rii + 0(77““^), when 77 ^ 0 and 77 ^ 00. However, one can 

similarly adapt to equations of state with quite general asymptotic behavior: assume that 77(77) behaves 
like 77(77) = 770 + vo{r]) and 77(77) = 777 + 777(77) (where vo{r]) and 777(77) can go to zero in different, arbitrary 
ways). Choose co = f{r]), where f{r]) can be any smooth, monotonically increasing function adapted to the 
asymptotics of 77(77), a suitable choice of f{r]) is one where f{r]) asymptotically satisfies /{rj) = 7/0(77)^ 
(fc < 1 ) and f{r]) = 777(77)“^. By this choice n{oj) becomes a differentiable function. Using the variables 
(77, q, oj) we obtain again where the duj/d^ equation now reads duj/d^ = —(7(77 1 ^)(/ ^(tl^))- Thus the 
system is similar to (i , and allows similar treatment. 


Static spherically symmetric perfect fluid models are also of interest in General Relativity. Relativistic 
stellar models with a variety of equations of state have been studied in some earlier papers, see, e.g, [ 0 , 
0 , |l 4 | . Several of the new ideas presented in this paper should also be useful in the general 
relativistic context, and might yield new results in that area. Note that in the low pressure limit 
the relativistic equations approximate the Newtonian equations, and thus the fixed point analysis of a 
relativistic dynamical system partly resemble the one presented in this paper, at least in the low pressure 
region (compare with |ll 4 | ). Moreover, it might be possible to derive relativistic inequalities by adapting 
the idea of using certain Newtonian polytropic solutions to divide state spaces into tractable regions. 


It is also possible to go beyond the perfect fluid assumption. In a future paper we will consider the 
Vlasov-Poisson equations for a self-gravitating collisionless gas. The corresponding dynamical systems 
formulation displays many similarities with the present situation, although new features also occur. 


Finally, it might be fruitful to go beyond the static and spherically symmetric assumptions and use the 
dynamical systems approach as a starting point for perturbation theory. 
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A Proof of theorems 

A.l Dynamical systems 

Let us briefly recall some main tools from the theory of dynamical systems (see e.g., Q, |Q). 

Definition, (u-limit set). Consider a dynamical system on R’". The u-limit set u^x) of a point x G R™ 
is defined as the set of all accumulation points of the future orbit of x. 

Remark. Correspondingly, the a-limit set a{x) is defined as the set of accumulation points of the past 
orbit. The w-limit of a set S' C R™ is w(S) = Ua; 6 S^(^)- The w-limit sets (a-limit sets) characterize the 
future (the past) asymptotic behavior of the dynamical system. The simplest examples for limit sets are 
fixed points and periodic orbits. An orbit whose lo- and a-limit set is a fixed point is called heteroclinic; 
a homoclinic orbit originates from and ends in one and the same fixed point. 

Theorem A.l. (LaSalle principle). Consider a dynamical system y = f{y) on R™ and the closed, 
bounded, future invariant set S C R™. If Z : S IS. is a function which is monotonic along orbits in 
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S, i.e., Z < 0 (Z >0), then 

w(S') C {s G S' I Z = 0} . (36) 

Replacing ’’future invariant” by ’’past invariant” results in the analogous statement for the a-limit set. 

Theorem A.2. (Monotonicity principle Consider a dynamical system on R."* and the invariant 

set S C R™. If Z : S ^ M. is a function which is strictly decreasing along orbits in S, then 

^{S) C {s G S\S I lim Z{x) yf sup.Z} (37a) 

x^s g 

a{S) C {s G S\S I lim Z{x) ^ inf Z} . (37b) 

x—>s S 

The monotonicity principle gives information about the global asymptotic behavior of the dynamical 
system. Locally in the neighborhood of a hxed point, the flow of the dynamical system is determined 
by the stability features of the fixed point. If the fixed point is hyperbolic, i.e., if the linearization of 
the system at the fixed point is a matrix possessing eigenvalues with non-vanishing real parts, then the 
Hartman-Grobman theorem applies: In a neighborhood of a hyperbolic fixed point the full nonlinear 
dynamical system and the linearized system are topologically equivalent. Non-hyperbolic fixed points 
are treated in center manifold theory. We give a brief introduction (see, e.g., or ||]): 

Consider an autonomous dynamical system y = F{y), where y{t) G R™. Assume that there exists a fixed 
point P and assume that the linearization of the system at P is described by a matrix which can be 
diagonalized. Accordingly, the phase space R™ can be decomposed into the direct sum if® © if“ © if^ of 
the stable subspace if®, the unstable subspace if“ and the center subspace E‘^. Moreover, without loss 
of generality, we can choose variables such that the dynamical system can be written as 


xi = AixiNi{xi,X 2 ) (38a) 

X2 = A 2 X 2 + N2{xi,X2) ■ (38b) 

Here, xi G if'’ is an 77ic"dimensional vector, 0:2 G if® © if“ is (mu + ms)-dimensional, the Ai are matrices 
representing the linearized system, and the Ni denote the nonlinear terms. 

A center manifold M” is an invariant manifold in the phase space R"*, i.e., a submanifold of R''", which is 
invariant under the flow (ft of the dynamical system. Moreover, M" must contain P and is required to be 
tangent to if". In a neighborhood of P one can describe M" as the graph of a function h : E^ ^ E’’(BE'^, 
i.e., {xi,h{xi)) G M'’. The function h is a solution to the differential equations 

dj;^h{xi) [Aixi + Ni{xi,h{xi))] = A 2 h{xi) + N 2 {xi,h{xi)), (39) 

which also satisfies h(0) = 0 and (i9a;j/i)(0) = 0 (tangency conditions). Recall that P = (0,0) in the 
adapted variables we use. The Shoshitaishvili theorem (reduction theorem, i, 0) generalizes the 
Hartman-Grobman theorem for the present situation: 

Theorem A.3. (Shoshitaishvili theorem). Consider the dynamical system in a neighborhood of the 
fixed point P = (0, 0). Then the flow of the full nonlinear system and the flow of the reduced system 

ii = Ai^i + Ai(a;i, h(a:i)) (40a) 

±2 = A 2 X 2 ■ (40b) 

are locally equivalent, i.e., there exists a local homeomorphism on phase space, such that = 
o ^reduced ^ Here, h is given by (0). 

Equation ( 0) is in general a nonlinear partial differential equation and cannot be solved in closed form. 
However, it holds that if a function h solves (pm up to 0(a:") as well as the tangency conditions, then 
h{xi) = h{xi) + 0{xf ). Hence, one can solve (|39D approximately by a formal power series. 
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A.2 Attractors in the state space 


In this subsection we prove that the global asymptotic behavior of the dynamical system ( 0 ) is fully 
determined by the fixed points of Table 0 and the closed orbits Ci, (72 (see Section 


Theorem |4.1 

to 5. 

X oo(—oo). 


All solutions converge to fixed points when X ±oo except when no (ni) is equal 
In this latter case solutions also converge to the 1-parameter set of closed orbits Ci (C 2 ) when 


Proof. All orbits of the dynamical system (Q) are strictly monotonically decreasing in O except when 
n = 0, n = 1, (3 = 0, and (7 = 1. It follows from the monotonicity principle that the w-limit set and the 
a-limit set of the state space reside on these invariant subspaces. 

More precisely, if x is an interior point of the state space, then 

a{x) C {P = 1} U {Q — 0} and uj^x) C {P = 0} . (41) 


To see this, note that (7 is a monotonically decreasing function in a large part of the state space including 
{3/4 < (7 < 1} X [0,1]^. This excludes attracting sets in this region and prevents interior orbits of having 
w-limits on the {U =1} side face. Considering the known orbits and the local structure of the fixed 
points on {Q — 0}, and making use of the Shoshitaishvili theorem AA, it becomes clear that the fixed 
points act as a-limit sets only, and not as w-limits. Finally note that Q is a, monotonically increasing 
function when 0 < (3 < 1 and 1/2 < [/, ensuring in particular that L4 cannot contain an a-limit of any 
interior orbit and thus leaving {fl = 1} U {Q = 0} as the only set containing a-limits. Note that if x lies 
on any of the side faces, then a{x),u!{x) are fixed points on the side faces themselves, as follows from the 
known solution structure. 

We now come to the remaining point and investigate what type of attractors can be found in the {fl = 0} 
plane (and the identical {fl = 1} plane). For uq < S the function 


Z = 


U 


1 -Uj \l-Q 


Q 


(42) 


is monotonic, which is seen from 


— = (2[/(l - Q) + (3 - n{n))Q{l -U))Z . 

CLA 


(43) 


This excludes attractors in the interior of the = 0 set when uq < 3. Consequently, combined with the 
solution structure on the boundary it follows that the set of possible attractors is represented by the fixed 
points discussed in Section ^ (it is easy to see that the equations prevent the existence of a heteroclinic 
cycle on the boundary when 0 < uq < 3). 

In the case np > 3 (ni > 3) we notice the fixed point P 3 {Pq) in the interior of the plane. To prove the 
fact that this is the only interior limit set as long as np 5 (ui 5) we resort to the simpler intermediate 
dynamical system (^ in the non-compactified variables {u,q). Define 


(k = (ug)"o-i 


no - 5 


-k 


2q 


{no - 1)^ np - 1 2 


qu X 
np -k 1/ 


(44) 


Note that for the special case np = 5, apart from the different choice of variables, (k -k (7 = 0 coincides 
with $ = 0 (see (^)), i.e., for np = 5, 'k -k (7 = 0 describes the orbits in the polytropic subset. 

More generally, for all np > 3 the equation (k -k (7 = 0 with (7min(np) < (7 < 0 represents a family of 
closed curves with the fixed point P3 as the center, which itself is represented by (7 = (7min. The shape 
of the (7 = 0 curve depends on np. For 3 < np < 5 it originates at some point below Pi on the (3-axis, 
turns around P3 and ends in a point above Pi on the (3-axis again. Increasing ng to np = 5 this (7 = 0 
curve becomes the regular orbit P2-P1 (plus the sections of the Q- and U-ax.es closing this curve, creating 
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^2 L 4 Pi Li Pi Li 



(a) no = -1 


(b) no = 5 


(c) no = 8 


Figure 12: The curves 'I' + C = 0. The polytropic indices no are: (a) uq = 4, (b) no = 5, (c) no = 8. 


a heteroclinic cycle). For hq > 5 the (7 = 0 curve originates from Li and ends in some point below Pi on 


the (5-axis. The remaining curves (7 > 0 always cover the second half of the = 0 plane. See Figure 12 


The function 4' is monotonic on orbits for no ^ 5 (ni ^ 5), i.e., 


— = (no- 5) ((no - 1) ^(uq) "o-i (-2 -h (no - l)g)^) . (45) 

Using the LaSalle principle, this implies that the limit-sets lie on g = 2/(no — 1). Taking, e.g., dq/d^ into 
consideration, this guarantees that P 3 (Pk) is the only possible tu-limit set (a-limit set) in the interior of 
the {U = 0} {{n = 1}) plane. 


Together with the solution structure on the boundary, the above arguments also imply that there are no 
heteroclinic cycles when no ^ 5 (ni ^ 5), and thus it follows that all orbits start and end at fixed points, 
with the exception of when no = 5 (ni = 5) when the closed orbits Ci (C 2 ) constitute a limit set. ■ 
Remark. The expression (Q) can be derived relatively easy in a set of variables different from the ones 
previously used, and then transformed to (u,q). Those variables are obtained from the original p, p, or 
77 by the so-called Emden transformations (see [||, p.90) 


A.3 The bifurcations at no = 3 and no = 5. 


Applying center manifold theory we prove that in continuation with no < 3, Pi is excluded as an 
attractor for interior orbits (and interior orbits on = 0) when no = 3. These considerations are 
intimately connected with the result that solutions have finite radii in the case no = 3 (see Theorem 3.1). 


To establish the claim it is preferable to investigate the dynamical system in the uncompactified variables 
(u, q,u!), i.e., the system (^. In this setting Pi is represented by the fixed point (0,1, 0); the eigenvalues 
of the linearization of the system at Pi are {3 — no,l,— a} and the corresponding eigenvectors are 
{(2 — no)ei -I- 62 , 62 , 63 } (where the form an orthonormal basis in the (n, g,a;)-space). Introducing 
adapted variables {x > 0 ,y, z > 0 } according to 


u = x(no — 2) q = 1 — X + y uj = z 


(46) 


we obtain the following dynamical system: 

X = (3 — no)x + a:[3 — n(z) — x(no — 2) — (y — x)n(z)] 

y= y + y'^ + x'^(5 + n{z) - 2no) - xy{4 + n(z) - no) - x(n(z) - no) (47) 

z = —az + a(x — y)z 
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Setting no = 3 we observe that ((1,0,0)) is the center subspace E‘^ of the system, ({(0,1,0), (0,0,1)}) 
corresponds to E^ (B E^^, so that our system ( ^ is of the form (|3^). 

Solving (^), as indicated in Subsection |A.l| , we see that the center manifold of the system ( ^ ) can be 
represented by the following local graph h{x) = {hy{x), hz{x)): 

hy{x) =—2x^ — 16x^ + 0{x'^) hz{x)=0. (48) 


By Shoshitaishvili’s theorem 
to the reduced system 


A.3 


in a neighborhood of the fixed point, the full system (^) is equivalent 

X = 2x^ + 6x^ + 0{x'^) 

y = y (49) 


z 


—az . 


Consider now a solution passing through a point {x > 0, y, z > 0) (corresponding to an orbit in the 
interior of the cube). Obviously, such a solution cannot converge to Pi = (0, 0,0). Accordingly, the claim 
is proved. 

It can also be shown that P 4 is the a-limit set for a two-parameter family of orbits in the case ni = 3. 
The investigation of P 4 is almost identical to the above one for Pi. 

We turn to the bifurcation at no = 5 and prove that each closed orbit Ci {C 2 ) acts as an w-limit (a-limit) 
set for a one-parameter family of interior orbits of the state space. The method of investigation is again 
based on center manifold theory. 

In a neighborhood of P3 (Pq) the full linear system is equivalent to the decoupled system 


i = fi{x,y) 

y = f2ix,y) (50) 

z = —az , 

where x,y are variables linearly depending on U,Q and z = O. Since we know from (^ that for z = 0 
the orbits are periodic we obtain the desired result. Note also that, since dn/dPl = 0, when fl ^ 0, the 
above type of decoupling also takes place asymptotically, when the orbits are not in a neighborhood of 
P3 (and analogously for Pg). 


A.4 


Proof of Theorems 6.4 and 


6.5 


Theorem |6.4| . (Spiral structure of the {M, R)-diagram). Consider an asymptotically polytropic equation 
of state with asymptotic indices 0 < no < 3 and ni. Then the mass-radius relation for high central 
pressures possesses a spiral structure, if and only if ni >5 with the spiral given by 


( R{Vc) 

V 


Ro 

Mo 



B J{'j2log — )b o{ 

Vc 



(51) 


where Rq and Mq are constants, B is a non-singular matrix, and b a non-zero vector. The matrix 
J{^p) G SO(2) describes a rotation by an angle (p, and the constants 71 and 72 are given by 


71 = (ni - 5) i 72 = i . 


(52) 


Recall from Table [J that Vb = \/—i — 22ni -|- 7nf. 

Sketch of proof. The proof is based on the idea that the spiral structure, which is formed by the orbits 
in a neighborhood of Pq , is reflected in a spiral in the (M, P)-diagram. 

Part 1: We show that ni > 5 implies a spiral in the (M, P)-diagram and conversely. 
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First assume that ni > 5, and note that in this case there exists one single orbit originating from Pq. 
Like every solution for an equation of state with no < 3, it possesses finite mass, Mo, and finite radius, 

Ro- 

Set rie = 1 — e, where e is chosen sufficiently small, so that n{Vl) > 5 for all Consider the = fig 

plane, and let I be the intersection point of the Pe-orbit with this plane. Choose a small neighborhood 
Uj of I within the plane. If e and Uj are small enough then, by the approximate decoupling of the 
dynamical system, the intersection of the one-parameter set of regular solutions, possessing sufficiently 
high central pressures, forms a spiral in [//. 

As for equation consider a small cylinder centered around P 2 - The Pe-orbit intersects the cylinder 
in a point O, with coordinates (ho,(t)o)i determining the mass Mq and the radius Rq- By the flow 
of the dynamical system we transport Ui to a neighborhood Uq of O on the cylinder. Provided that 
Uj has been chosen sufficiently small, the spiral in Uj is mapped to a spiral in Pq, since the map is a 
diffeomorphism. 

Inserting h = ho + Sh and (j) = (j)o + Scj) into (p0|), leads to the following mass and radius relations in Uq', 


R = Ro + ciSh + c 2 S(j) M = Mo + coSh + , (53) 

where the vectors (ci, C 2 ), (ca, C 4 ) are linearly independent. Accordingly, a spiral in results in a 

spiral in the (M, P)-diagram. 

Conversely, assume that the high pressure regular solutions form a spiral structure in the (M, P)-diagram. 
We then need to show that ni must be greater than 5. 


To this end, let {Mo,Ro) be the center of the spiral and O = {ho,(t>o) be the associated point on the 
cylinder near P 2 . Choose a small neighborhood Uo of O on the cylinder. By following the orbit passing 
through O backwards in A it can be seen that the (indirect) assumption ni < 5 leads to a contradiction. 
Namely, there are three possible origins of the considered orbit: If the O-orbit comes from Pg (plus C 2 in 
the case ni = 5), then this must be true for all orbits passing through Uo by continuous dependence on 
initial data. Accordingly, none of these solutions is regular, which is a contradiction to the assumption 
that the regular solutions form a spiral in Uo ■ The same argument holds if the O-orbit originates from 
Li. Finally assume that the orbit corresponds to a regular solution that comes from Oc on L 2 - However, 
not all regular solutions Q > Qc can pass through Uo, once more by continuous dependence on initial 
data, which is again a contradiction to the assumed behavior of the high pressure solutions. Thus the 
first part of the theorem is established. 

Part 2: We now give a sketch of the proof for the qualitative form of the spiral structure (|3ll). 


Consider again the orbit that originates from Pg and passes through the points I G Uj C {O = and 
O G Uo Q cylinder. For initial data lying in Ui close to I we may linearize the full dynamical system 
along the Pg-orbit; we get 

/ SU{X) \ / SUo\ 

SQiX) = A(A) ^Qo . (54) 

V 6n{x) J V 0 ; 


Below we use equation (54) to map the spiral of regular solutions {SUo{0,c), SQo{0,c)) (parametrized by 
the central value lie) down towards the cylinder around P 2 , and thus obtain the (M, P)-spiral. The next 
step is therefore to derive an expression for the spiral of regular solutions in Uj. 


Let Up{X), Qp{X) denote the regular solution of the polytropic system (^-|-(7b) with polytropic index 
ni- This solution gives rise to a spiral whose asymptotic form in the limit X 00 reads 


{ Up{X) 

V Qp(A) 


Up, 

QPe 


+ exp(-(5iA) B'Ji 62 X) b', 


(55) 


where Up,,Qp, are the coordinates of the fixed point Pg; B' is a non-singular matrix, and h' a non-zero 
vector; i5i = f (u-i — 5) and 82 = jVb. The constants Si and S 2 are the real and imaginary part of the 
complex eigenvalue that is associated with Pg, compare with Table |^. 
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Now consider the dynamical system for SI G Since we can choose e arbitrarily small (and 

since in addition = 0) we are justified to approximate n(Sf) by n(Sd) = ni in so that the 

system decouples. Accordingly, can be reduced to = «(! — S^)Qp(A)(l — t/p(A)), which can 

be solved: 

1 - Sf(A) = (1 - Oc) exp J Qp(A)(l - {7p(A))dAj , (56) 

where Sic is the initial value on L 2 . Inserting ( |5^ into (|^) we obtain the approximation 

1 — S1(A) = (1 — Sic) const exp (^aQpg{l — Up^) {X — X)), (57) 

where A and const are independent of Sic. 

We conclude that a regular solution that starts from Sic near 1 reaches [// C {SI = Slg} at Ac(Slc) = 
log yrn—^ const, where const is independent of Sic. Note that f3 = Qpf^{l — Up^) (compare with 
Table |^) . Transforming from SI to 77 yields 


Ac(??c) = ^ log?7c + const. 


(58) 


At A = Ac the ([/, (3)-components of the regular solution are located at some position on the spiral (55), 
so that we are able to explicitly give the spiral of regular solutions in [7/ C {SI = Sic} by inserting 
into (1^ : 

SUoiVc) \ ^ (^Y 
SQaiVc) j KVc) 

where B" and b" have absorbed the constants and a fixed rotation matrix. 


^(72 log 77 c ) 6" 


(59) 


Inserting ( p^ ) into (p^ and projecting onto the small cylinder at P 2 affects the spiral by an additional 
distortion as does the change of coordinates from {h, (p) to (M, R). Hence Theorem ^ is established. ■ 

Remark. Theorem |6.4| can easily be generalized to the case no = 0; however, as all orbits end in L 5 instead 
of P 2 , the (M, i?)-relations dil) and ( po|) must be replaced by different ones. Clearly these relations are 
not qualitatively different. In order to visualize things, one can again make use of the cylinder picture: 
consider a cylinder lying down with small radius with L 5 as its center. Then every orbit can be uniquely 
characterized by its intersection with this cylinder, i.e., by ([/o,</>). Accordingly, all arguments used in 
the proof of Theorem go through. 


Theorem 


6.5 


(Polytropic behavior of the (M, R)-diagram). Consider an asymptotically polytropic 
equation of state with asymptotic indices 0 < no, and 0 < rii < 5. To first order, the mass-radius 
relation for high central pressures is approximated by the mass-radius relation for an exact polytrope with 
polytropic index ni, i.e., 


R{Vc) = Rp-.m iVc) C[n{v)] M ( 77 ^) = ( 77 ^) C[n{r])] 


(60) 


where R^.,ni(r]c) and are the radius and the mass for an exact polytrope with index ni (compare 

with i.e., 


Rp:ni(r/c) = Cii(7li)?7, 


1—ni 


A7p:„i(77c) = CMini)r], 


3 — 711 


and C is a functional acting on the index-function nfq), given by 

C[n{r])] = exp I i J 


1 00 

n{r]) - no , 1 f 77(77) - 77i 

- dV+o - 

V 2 7 77 

0 1 


dr] 


(61) 


(62) 
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Sketch of proof. The proof of the claim is quite analogous to the proof of Theorem For high central 
pressures (i.e., close to = 1 } in the state space) the dynamical system almost decouples (recall 
that dn{0)/dPl = 0 at = 1). Therefore, for A smaller than some value A, the (U, Q) component of orbits 
that originate from close to 1 , can be approximated by the corresponding regular polytropic 

orbit (C/p, Qp) in the {fl = 1} plane. Consequently the O component behaves like 


f2(A) = 1 — (1 — flc) exp 



Qp(A)(l-C/p)(A))dA 


(63) 


Choosing Q.c sufficiently close to one, the corresponding orbit will come arbitrarily close to P 5 at A = A. 
Thus, linearizing the dynamical system along the P 5 -P 2 orbit, we can calculate the further evolution of 
the high pressure orbits. We obtain 


C7(A) = Pp(A)e-^(’'(^»e^(''(^» 

(64a) 

Q(A) = l-(l-Qp(A))e-‘ 

(64b) 

0 {\) = (! + (!- r!(A))e-“^e“^)"^ , 

where N{r]) is the integral 

■=- , 

J V 

(64c) 

(65) 

and 77 (A) corresponds to n(A), i.e., 77(A) = n(A)^/“e^e“'^. 


For large A, when n(A) (respectively 77 (A)) has become small, we can 
pressure solutions by 

therefore approximate the high 


C/(A) = Pp(A)e"«^e-"«^ exp ^ 

Q(A) = l-(l-Qp(A))e-* ( 66 b) 

0(A) = (1 - 0(A))-ie“^e-“^ , ( 66 c) 

where we note that the last exponential in (| 66 a|) can be simplified by using 


no 


dr] 


( 66 a) 


n{T]) - no 


dr] = (m - no)(log? 7 c - / (5p(A)(l - Pp)(A))dA) x 

J —00 


n{r]) - no 


dr] 


i{r]) - ni 


dr]) . 


(67) 


Equations ( |^ exhibit the typical form of linearized solutions near the hxed point P 2 . Recall from ([l9| ) 
that we can easily read off the total mass and radius of such a solution, yielding 

R^ (X Tyi-"! oc . (68) 

Moreover, comparing constants appearing in the general expressions for M and R and the analogue 
expressions for exact polytropes, we see that they differ merely by a factor expressible as C[?^(? 7 )], which 
completes the proof of the theorem. ■ 
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B Miscellaneous generalizations 


B.l The negative mass cube 


Also perfect fluid solutions with negative mass, i.e., solutions of (^) with negative m(r) can be studied 
in our setting. As before we use compactified variables {U,Q,Ll) 


U = 


u 

TTR ’ 


Q 


q 

i + M ’ 


Ll = 


UJ 

1 + Cl; ’ 


(69) 


where (u, q, oj) are defined as in (Q). By our assumption m < 0 it follows that LI is positive and that U, Q 
are negative. Including the boundaries in our discussion we obtain the state space [—1, 0]^ x [0,1]. 

Introducing the independent variable A according to dX/d^ = (I + C/)“^(I + <5)“^ yields the following 
equations for U, Q, LI: 


— =U{1 + U)[{1 + Q)(3 + 2U) - n(Ll) Q{1 + U)] 
ctA 

= Q(1 + Q)[—(I + Q) + Q(1 + U)] 

^ = -an{i-n)Q{i + u ). 


(70a) 

(70b) 

(70c) 


Compared to the positive mass cube the negative mass cube possesses a much simpler fixed point structure 
(see Table ^) . 


Fixed point 

u 

Q 

Ll 

Eigenvalues 

Restrictions 

Li 

-1 

0 

LIq 

-5 , -1 , 0 


L 3 = L 3 

0 

0 

LIq 

3 , -1 , 0 


u 

-1 

-I 

LIq 

0 , 0,0 


L 5 

Uo 

-I 

0 

0 , (1 + Uo) , a(l + Uo) 

0 

II 

0 

^6 

Uo 

-I 

1 

0 , (1 + Uo) , -a(l + Uo) 

7 Zl = 0 

A 

0 

-I 

0 

TT-O 7 1 ; ^ 


h 

0 

-I 

1 

ni , 1 , —a 



Table 4: Local properties of the fixed points of the negative mass cube. 

For orbits in the interior of the state space the flow of the dynamical system is characterized by a 
universal a-limit set, namely the fixed point P 2 and a universal w-limit set, namely Li (in the special 
case no = 0, P 2 must be replaced by L 5 ). In brief: all orbits originate from P 2 and end in Li. 

Accordingly, using the linearization of ( [^ around P 2 and Li, we can easily derive the qualitative features 
of a negative mass solution: In a neighborhood of {r < i?_} the negative mass solution is described by 
a vacuum solution with the potential v = —M-jr + + vs, where AM is negative and V- = M-/R-. 

At r = i?_ the quantity 77 (which was defined as 77 = —v + vs) and therefore LI has its zero; here the 
solution enters the interior of the negative mass cube (at P 2 ), 

v{r) =—CSr + vs + 0{Sr^), (71) 


where Sr = and C = 

For r > i?_ the density is positive so that the solution begins to acquire ’’positive mass”. This happens 
continuously until the radius r = R is reached where the mass has been reduced to zero, m{R) = 0. At 
r = R the potential v(r) has the (negative) minimum vq + vs- In the state space this corresponds to the 
orbit reaching Li at a certain LIq- 
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The solution can be continued into the positive mass cube by identifying Li with Li: starting from S7o 
on Li the solution proceeds as a solution which now possesses positive mass. 

Equation ( 0 ) gives the approximate form of the solution near the radius R. Note that this expansion 
holds for both Jr < 0 and Jr > 0. 

B.2 Non-increasing eqnations of state 

In Section ^ we require no to be a non-negative constant (where the special case np = 0 corresponds 
to an equation of state with p const > 0 when p ^ 0). (The case no < 0 could be interpreted as 
p —t oo when p ^ 0). However, the non-negativity requirement is not necessary for the index-function 
n{r]) (or n(Ll)) as a whole, however, recall that a non-negative index-function n(ry) is associated with a 
monotonically increasing equation of state p(p). 

In particular, ni is allowed to be negative, so that in addition to the patterns of Figure]^ the orbits on 
the {fl = 1} plane can look like in Figure 


Pi Li 



Figure 13: Orbits on the H = 1 plane in terms of the variables U and Q for negative values of ni. 

Note that although Figure qualitatively differs from the polytropic pictures (Figure ^), the global 
asymptotics for interior orbits does not change. 

In order to prove that a closer investigation of L 4 is necessary, namely, we need to show that no interior 
orbits converges to L 4 for A ^ ± 00 . This can be done using polar coordinates (r, (p) centered at L 4 and 
a changed independent variable. 

B.3 Superposition of solutions 

In this subsection we investigate superpositions of perfect fluid solutions. 

Recall that rj{r) = vs — v{r), where v is the Newtonian potential. The equations of the self-gravitating 
perfect fluid are the Euler-Poisson equations leading to 

Ary = —drrp , (72) 

where p is to be understood as p{r]), the equation of state. This formulation is equivalent to the system (^. 

If rii{r) (r = 1, 2) solves 0 ), Arji = —Airpi, then the superposition 77 = ryi -f 772 -f c (c a constant adapted 
to the boundary conditions) satisfies Ary = — 47 r(pi( 77 i) + p 2 (? 72 ))- In this case there exists an equation 
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of state such that Arj = —4:Trp{fi). Clearly, this equation of state p{fi) depends on the explicit form 
of Pi (r) and p 2 (r). In this (weak) sense the superposition of two perfect fluid solutions is again a perfect 
fluid solution. 


A particularly interesting case to study is the superposition of a regular perfect fluid solution with a 
vacuum solution corresponding to a point mass: 

Consider an asymptotically polytropic equation of state p{rf) with asymptotic behavior given by the 
indices uq and ni and consider an associated regular perfect fluid solution r]{r). Let us assume that it is 
of finite extent with radius R and total mass M. A point mass solution is given by the potential —Mq/t, 
where Mq is some positive constant. It satisfies the vacuum equation on {r > 0}, i.e., A{MQ/r) = 0. The 
superposition of these two solutions is given by fj{r) := ri(r) + Mq/t — Mq/R. fjir) satisfies the equation 
Afj = —Lkp(j]) and the boundary condition f)(i?) = rj{R) = 0. 

Note that pir) is monotonically decreasing; hence the same is true for f]{r) and there exists the inverse 
function r{f)). Define 

/i«) :=/.(.)) = />(«-4| + , (73) 

then Afj = —I'Kp{f}) (on {r > 0}), i.e., fj{r) is a perfect fluid solution corresponding to the equation of 
state pip). By using the expansions 0) and of p{r) it can easily be shown that the equation of 
state p(ri) is asymptotically polytropic with asymptotic indices ho = no and hi = 0, namely 

Kv) ^ P- \ m + Mo ) + (^^0) > PiR) = Pc{l + Oirj-^)) {p ^ oo) . (74) 

In particular we see that pifi) has the same asymptotics as p{p) for low densities and behaves like an 
incompressible fluid for large pressures. 


In the dynamical systems formulation the perfect fluid solution p(r) appears as a solution that is emitted 
from the fixed point P4 in the state space associated with p{p) (compare with (12a) for ni =0). We can 
take the viewpoint that, by superposing a regular solution p{r) (associated with p{p)) with a mass point 
solution, the regular solution is mapped to a solution emitted from P 4 (in the p(r7)-state space). 


Remark. As an example we have taken regular solutions with finite extent. However, the statement is the 
same for general regular solutions. Moreover, maps as described above can be obtained for completely 
different families of solutions in a similar way. 


B.4 Modified formulation 


As we have already stressed earlier, the formulation presented in Section in particular the choice of 
variables (Q), is adapted to asymptotically polytropic equations of state in the proper sense, i.e., n < 00 . 
However, in some situations one is also interested in equations of state p{p) that behave asymptotically 
like p (X p'^°, p oc p'^^ for arbitrary vo,vi > 0. In particular one would like to cover equations of state 
that are asymptotically linear in the low and high pressure regime (i.e., r'o = 1 or = 1). 

In order to deal with this slightly more general class of equations of state, we can introduce the variables 


u = 


Anr^ p 


V = 


mp 


a; = 


m rp 

and the corresponding compactified C/, V, LI. In these variables we obtain the dynamical system 


— = C/(I - [/)[(! - H)(3 - 4C/) - n{n) H(I - U)] 

^ = 1/(1 - V)[{2U - 1)(1 - H) + (1 - nm y(l - U)] 
^ = -oH(l-H)y(l-C/), 


(75) 


(76a) 

(76b) 

(76c) 
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where ^{O,) = iy{p{Ll)) is the index-function defined via i'{p) = For polytropic equations of state 

this is a constant, v{p) = v = The index-function ^{p) is related to the so-called adiabatic index 

7 ad(p) via v(p) = l/ 7 ad(p) (see, e.g., |^). Note that for general equations of state v{p) ^ ri^p)+i (where 
n{p) = n{r]{p)) is the index-function defined in (^). 

The flow of this dynamical system on the = 0 plane is depicted in Figure 0 this is the counterpart of 
Figure in the present formulation. 



L3 L2 Li 


(a) = 0 


P2 L, 





P2 Li 



P2 Li 



(c) = I 


P 2 Li 



(g) = 1 


P2 Li 



(h) u = 1.5 


Figure 14: Orbits in the polytropic subset = 0 in terms of the alternative variables U (horizontal axis) 
and V (vertical axis). 

As indicated by Figure ^ all solutions have infinite radius and infinite mass for izq > 1 . This is in 
accord with the criterion guaranteeing infinite extent if 77 = dp'p~^{p') 00 (see for a good 
presentation). For further statements, the theorems discussed in Section can be reformulated and 
adopted to the present context. 

Compared to the formulation presented In Section |^, the formulation (|75| ) has certain disadvantages: 
Firstly, v appears twice in the dynamical system (0), secondly, the fixed points on the boundaries of the 
state space are not “fixed”, i.e., their position depends on Vi. Moreover, Chandrasekhar’s equation of 
state is quite complicated in terms of p,p variables, while it is quite elegantly expressed in p, rj variables. 
Thus we find that ’’ 77 -based” variables, in many contexts, is a quite natural choice. 
























REFERENCES 


36 


References 

[1] H. Buchdahl. General Relativistic Fluid Spheres. Phys. Rev., 116 : 1027-1034, 1959. 

[2] J. Carr. Applications of center manifold theory. Springer Verlag, New York, 1981. 

[3] S. Chandrasekhar. An introduction to the study of stellar structure. University of Chicago Press, 
Chicago, 1939. 

[4] J.D. Crawford. Introduction to bifurcation theory. Rev. Mod. Phys., 63-4 : 991-1038, 1991. 

[51 J.M. Heinzle. (In)finiteness of static spherically symmetric perfect fluids. Class. Quantum Grav., 
19 : 2835-2853, 2002. 

[6] G.P. Horedt. Topology of the Lane-Emden-equation. Astron. Astrophys., 177 : 117-130, 1987. 

[7] H. Kimura. A Study of Simple Polytropes I. Fundamentals and Classification of Solutions. Publ. 
Astron. Soc. .Japan, 33 : 273-298, 1981. 

[8] R. Kippenhahn and A. Weigert. Stellar structure and evolution. Springer, Berlin Heidelberg, 1994. 

[9] L. Lindblom. On the symmetries of equilibrium star models. In S. Chandrasekhar, editor. Classical 
General Relativity, Oxford Science Publications, pages 17-28. Oxford University Press, 1993. 

[10] T. Makino. On spherically symmetric stellar models in general relativity. J. Math. Kyoto Univ., 
38-1 : 55-69, 1998. 

[11] T. Makino. On the spiral structure of the (R,M)-diagram for a stellar model of the Tolman- 
Oppenheimer-Volkoff equation. Funkcialaj Ekvacioj, 43 : 471-489, 2000. 

[12] T. Makino. Ordinary Differential Equations of Emden Type Funkcialaj Ekvacioj, 27 : 319-329, 
1984. 

[13] U. Nilsson and C. Uggla. General relativistic stars: Linear equations of state. Ann. Phys., 286 : 278- 
291, 2000 . 

[14] U. Nilsson and C. Uggla. General relativistic stars: Polytropic equation of state. Ann. Phys., 
286 : 292-319, 2000. 

[15] A.D. Rendall and K.P. Tod. Dynamics of spatially homogeneous solutions of the Einstein-Vlasov 
equations which are locally rotationally symmetric. Glass. Quantum Grav., 16 : 1705-1726, 1999. 

[16] A.D. Rendall and Bernd G. Schmidt. Existence and properties of spherically symmetric static fluid 
bodies with a given equation of state. Class. Quantum Grav., 8 : 985-1000, 1991. 

[17] U.M. Schaudt. On Static Stars in Newtonian Gravity and Lane-Emden Type Equations. Ann. Henri 
Poincare, 5 : 945-976, 2000. 

[18] S.L. Shapiro and S.A. Teukolsky. Black Holes, White Dwarfs and Neutron Stars. John Wiley & 
Sons, 1983. 

[19] W. Simon. Criteria for (in)finite extent of static perfect fluids. In J. Frauendiener and H. Friedrich, 
editors. The conformal structure of space-times: Geometry, Analysis, Numerics, Lecture Notes in 
Physics 604, Springer (2002). 

[20] J. Wainwright and G.F.R. Ellis. Dynamical systems in cosmology. Cambridge University Press, 
Cambridge, 1997. 



